Signature of a Fano-resonance in a plasmonic meta-molecule's local density of optical states 
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We present measurements on plasmonic meta-molecules under local excitation using cathodoluminescence 
which show a spatial redistribution of the local density of optical states (LDOS) at the same frequency where a 
sharp spectral Fano-feature in the extinction cross section has been observed. Our analytical model shows that 
both near- and far-field effects arise due to interference of the same two eigenmodes of the system. We present 
quantitative insights both in a bare state, and in a dressed state picture that describe plasmonic Fano interference 
either as near-field amplitude transfer between three coupled bare states, or as interference of two uncoupled 
eigenmodes in the far field. We identify the same eigenmode causing a dip in extinction to strongly enhance the 
radiative LDOS, making it a promising candidate for spontaneous emission control. 



Introduction. Interference is ubiquitous in physics. Sig- 
nificant recent advances in optics as well as quantum physics 
hinge on interference, inherent in the wave nature of light 
and matter, and the superposition principle. In quantum op- 
tics, the Fano-effect and its occurrence in electromagneti- 
cally induced transparency (FIT) have in particular triggered 
tremendous interest as phenomena relying on quantum inter- 
ference Q in light-matter coupling. In FIT, a strongly absorb- 
ing atomic vapor coupled to an intense pump field acquires a 
narrow transparency window, with unusual features, such as 
ultralow group velocities and huge nonlinearities. These ex- 
traordinary properties have attracted the interest of the field 
of nano-photonics, the science of engineering the genera- 
tion, the propagation and the absorption of light on a sub- 
wavelength scale |2|. The aspiration of nano-optical circuitry 
with powerful functionality led to the development of opti- 
cal meta-materials. These artificial materials are composed 
of meta-atoms, designed building blocks giving rise to pe- 
culiar properties not found in natural materials |3]. Inspired 
by quantum optics, scientists have identified plasmonic meta- 
molecules whose optical properties mimic FIT-lineshapes 
in atomic vapors, an effect termed 'plasmon-induced trans- 
parency' (PIT) |4 |, based on the Fano-interference of a super- 
and a sub-radiant mode. Fven without the benefit of a full 
electrodynamic model reaching beyond brute force numerical 
simulations, remarkable intuition and simple electrostatic ar- 
guments have led to the development of several structures ex- 
hibiting PIT | 4 - 13 |. While in PIT plasmonic meta-molecules 
control the propagation of light by creating narrow dark reso- 
nances useful for slow light or sensing, another class of nano- 
structures termed 'optical antennas' is currently being devel- 
oped to tailor light matter interaction lfT4tiT6i . Antennas ex- 
ploit bright resonances to enhance the emission of light. Prac- 
tically all aspects of spontaneous emission control by optical 
antennas rely on designed enhancement of the local density of 
optical states (LDOS), arguably the most fundamental quan- 
tity in nano-optics |2|. An outstanding question is what the 
LDOS of PIT structures is, and if one can use the narrow dark 
lines of PIT to improve optical antennas. To answer this ques- 
tion it is essential to unravel which modes are involved in PIT, 
how they project on localized driving, and how they give rise 
to our observable, i.e. far field radiation. 




Figure 1. (a) SEM micrograph of plasmonic heptamer. Scale bar 
is 200 nm. (b) Sketch of experimental setup showing sample, im- 
pinging electron beam, parabolic mirror and optics guiding CL to a 
spectrometer. 

In this Letter, we map the LDOS of a plasmonic molecule 
known to exhibit PIT. Under near-field driving using cathodo- 
luminescence (CL) we observe a marked spatial redistribu- 
tion of LDOS occurring at wavelengths coincident with the 
PIT-dip in extinction, thereby relating far-field data to the 
near-field LDOS. We present a model that identifies the di- 
agonal 'dressed states' of the meta-molecule and proves the 
LDOS feature to be due to interference of the same eigen- 
states that are responsible for PIT in plane wave extinction. 
We draw three remarkable conclusions beyond the analogy 
with atomic FIT. Firstly, we find a 'screening' state that signif- 
icantly renormalizes the interaction. Secondly, we conclude 
that all observations can be viewed purely as far field inter- 
ference of two eigen-modes of the system. Thirdly and sur- 
prisingly, we identify the sub-radiant mode to be a promising 
candidate to enhance brightness and rate of spontaneous emis- 
sion into the far field. 

Experimental We fabricated plasmonic heptamers on a 
Si wafer covered with 1 jam thermal oxide by electron beam 
lithography |17|, thermal evaporation of 35 nm Au, and lift 
off. Fach heptamer consists of nominally identical particles 
arranged on the corners and center of a hexagon. A SFM im- 
age of a typical heptamer with particle diameter 150nm, and 
gap width 25 nm is shown in Fig.[TJa). The CL measurements 
were performed in a scanning electron microscope, sketched 
in Fig. [TJb). CL maps are acquired by raster scanning the 
focused electron beam (30 keV, waist < 5 nm) across the sam- 
ple. Fmitted light collected with a parabolic mirror (accep- 
tance angle 4.6 sr) is spectrally analyzed on a silicon CCD 
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Figure 2. (a) Measured CL excitability maps of plasmonic heptamer. 
(b) Calculated maps of radiative LDOS. Insert: colorbar for (a) and 
(b). (c) Limits of colorscale for (a). Background around 700 nm 
stems from oxide substrate, (d) Calculated radiative LDOS enhance- 
ment on central (purple) and outer (olive) particles of heptamer with 
particle diameter/gapwidth 150/25 nm. Blue and red symbols de- 
note colorscale limits of frames in (b). (e) Measured differential ex- 
citability spectra showing intensity difference between outer and cen- 
tral particles. Green (black): Particle diameter/gapwidth 150/25 nm 
(100/25 nm). Errorbars denote standard deviation of 4 (2) structures, 
(f) Calculated power radiated by super- and subradiant eigenmodes 
of heptamer upon driving central (purple)/outer (olive)/all (orange 
dashed) particles with vertical electric field. 



camera The signal hence represents 'CL excitability' as 
a function of detection wavelength and spatial excitation coor- 
dinate. In Fig.|2ja) we show CL collected from the heptamer 
in Fig. [TJa) as a series of spatial excitation maps by binning 
the data into 50 nm wavelength slices. We identified pixels 
on substrate vs. on Au particles by thresholding the SEM data 
collected in parallel with the CL. We clamp the color value for 
all substrate pixels to the smallest value obtained on any Au 
particle to maximize color contrast for the regions of interest 
we analyze, i.e. the particles. The minimum and maximum 
values of the colormaps in Fig. [2ja,b) [colorbar in insert in 
Fig. [2jb)] are plotted in Fig. |2jc,d), respectively. While at 
wavelengths shorter than 800 nm in Fig.|2ja) the outer parti- 
cles appear brighter, i.e. their excitability is higher, the sit- 
uation is reversed at wavelengths longer than 800 nm, where 
the inner particle is more excitable. Local hotspots are due 
to surface roughness inherent to the fabrication. To quantify 
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Figure 3. (a) Bare state non-diagonal polarization basis, consist- 
ing of a super-radiant, sub-radiant and dark symmetry adapted lin- 
ear combination (SALC) of Eiu derived from group theory, (b) 
Dressed state basis of decoupled superradiant, subradiant and dark 
eigenmodes (small out-of-phase components not visualized), (c) 
Real and (d) imaginary parts of polarizabilities of the eigenmodes in 
10~^^ Cm^/V. (e) and (f) are radiation patterns of superradiant and 
subradiant eigenmodes in W/sr. Radiative strength indicates total ra- 
diated power relative to superradiant eigenmode. 



the observed swap in local excitability we discriminate pixels 
belonging to outer or inner particles, spatially average their 
spectra and subtract them from each other. The resulting 'dif- 
ferential excitability' is therefore positive if the outer particles 
are more excitable than the inner one. This procedure elim- 
inates the broad background fluorescence of the oxide layer 
around 650 nm, which causes the overall increase in absolute 
signal towards the blue in Fig.[2jc). Note that this background 
together with the vanishing detector efficiency towards the in- 
frared cause the apparent discrepancy of Fig.|2jc) and (d). Fig- 
ure [2je) shows differential excitability spectra for heptamers 
with particle diameter/gap size 150/25(±5)nm averaged over 
four structures [green squares] and 100/25(±5)nm, averaged 
over two [black squares]. The wavelength where the differ- 
ential excitability changes sign, i.e. where excitability swaps 
from outer to inner particles, appears around 800 nm for the 
large heptamers, as reported in Fig. [2ja), and blue- shifts to 
650 nm for smaller particles. 

Electrodynamic model. We model the plasmonic hep- 
tamers with a fully electrodynamic coupled dipole model l|T9l . 
First, we show that this model captures both the far-field Fano 
interference, and the excitability swap in our experiment. Sub- 
sequently, we derive new insights by identifying the bare and 
dressed states in analogy to the atomic system. In our model, 
the response of each particle to an electric field, i.e. its dipole 
moment p = aE, is proportional to its Lorentzian polar- 
izability tensor a [20] . The resulting set of coupled linear 
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equations can be written in the form p = • Eg^ct, where 

F/ext is the incident field driving each particle. Both Eea^t and 
p have 3N=21 elements for N=l particles. The interaction 
matrix M has the inverse of a on the diagonal, while the off- 
diagonal elements describe the interaction between dipoles as 
set by the electrodynamic Green's function. For any dipole 
assembly M can be inverted to find the polarization state p in- 
duced by any driving Fiext • Since this linear problem contains 
full electrodynamic interactions, its solution allows to calcu- 
late any near- and far-field observable. Specifically, the inten- 
sity of generated CL is proportional to the component of the 
radiative LDOS along the impinging electron beam 12111221 . 

Figure [2|b) shows calculated maps of the radiative LDOS 
for driving perpendicular to and located 40 nm above the plane 
of the meta-molecule. As in the experimental data, the ring- 
shaped profiles per particle imply that we detect in-plane in- 
duced particle polarizations, despite the out-of-plane incident 
electron beam 1231 . The calculations show the measured 
swapping of LDOS from outer to inner particles with increas- 
ing wavelength [Fig|2jb)] as well as the shift of spectral fea- 
tures to shorter wavelengths with decreasing particle diame- 
ters [Fig. |2jf)]. As we show below, the Fano-dip in extinc- 
tion ifTTl calculated using the same model [see Fig. |4| coin- 
cides with the spectral position of the excitability swap in CL, 
pointing at a direct relation between LDOS and Fano interfer- 
ence. 

Symmetry of the heptamer and eigenmodes. While it is 
gratifying that our model quantitatively confirms both the 
Fano extinction dip and the concomitant CL signature, the 
analysis sofar provides little insight. Three essential steps 
simplify the problem. Firstly, as pointed out by Mirin et 
ah [24], the heptamer bears D6h symmetry. Therefore, Sym- 
metry Adapted Linear Combinations (SALCs) of dipole mo- 
ments 1251 cast the electrodynamic coupling matrix M into 
block diagonal form since group theory usefully extends be- 
yond electrostatic hybridization |26|, allowing a symmetry- 
based decoupling also for calculations of scattering. Secondly, 
since we rely on far field detection, we focus our attention 
on the infrared- active Ei^ irreducible representations. Due to 
the degeneracy of horizontal and vertical polarization in the 
sixfold symmetry we are left with a three-dimensional sub- 
space. Our choice of SALCs for Ei^ (vertical polarization 
only) are shown in Fig. [SJa). Two SALCs are the in and out 
of phase superpositions of one hexamer and the single particle 
Flu modes pointed out in 1241 . In addition to the superradiant 
S ALC (large net dipole moment) and subradiant S ALC (three- 
fold smaller dipole moment), symmetry requires a third dark 
SALC with zero net dipole moment. As a third essential step 
we note that while symmetry decouples different irreducible 
representations, we must still diagonalize the Ei^ submatrix 
of M to find the true decoupled eigenmodes. 

Remarkably, the three eigenvectors of Eiu, sketched in 
Fig. |3];b), are almost unchanged across the frequency range 
from 400 to 1000 nm. The complex eigenvalues show strong 
dispersion, as shown in Fig. [3jc,d). These eigenvalues are 
the 'eigenpolarizabilities' of the eigenmodes under which the 



Flu submatrix of M is diagonal and which therefore are by 
definition decoupled. Considering the dipole distributions in 
Fig. |3jb) we can classify the first eigenmode as superradiant, 
therefore featuring a broad eigenpolarizability [Fig.[3jc,d), red 
line], reminiscent of the superradiant SALC in that all dipole 
moments are approximately aligned [compare Fig. |3ja,b)]. 
The second eigenmode is subradiant with a narrow reso- 
nant eigenpolarizability around 630 nm, i.e. at the observed 
Fano dip and LDOS feature [Fig. [3jc,d), blue]. Due to the 
coupling set by IM, this mode is remarkably different from 
the sub-radiant SALC. The third mode has a narrow reso- 
nance at significantly shorter wavelengths, beyond our range 
of interest [green]. Since its eigenpolarizability is negligible 
around 630 nm, the Fano interference is explained in terms 
of just the super- and subradiant modes [red and blue in 
Fig.[3jc,d)]. These modes have largely overlapping radiation 
patterns [Fig. [3je,f)], although with a twelvefold ratio in in- 
tegrated radiated intensity at identical amplitude. While the 
eigenmodes are by definition decoupled, the excellent radia- 
tion pattern overlap implies that nearly completely destructive 
or constructive interference can occur in the far field. 

Two eigenmodes explain all observations. To demonstrate 
that two eigenmodes capture all the physics observed both un- 
der local and plane wave driving, we perform calculations us- 
ing just the super- and subradiant eigenmodes. The purple 
curve in Fig. [2jf) shows the radiated power when only the 
central particle is driven. The total radiated power (propor- 
tional to CL intensity generated on the central particle) shows 
an asymmetric peak at 650 nm with a steep slope on its blue 
side, and a strongly broadened wing on its red side. The local- 
ized driving mainly projects on the resonant subradiant mode. 
The superradiant mode provides a weak contribution that in- 
terferes destructively on the blue, and constructively on the 
red side of the resonance, leading to the typical asymmetric 
Fano-lineshape. Importantly, when driving the outer particle 
[olive curve in Fig.[2jf)], the peak asymmetry is reversed, as 
the relative phase between excitation of the broad superradi- 
ant and narrow subradiant mode is swapped. This asymmetric 
broadening in opposite directions underlies the measured dis- 
persive differential excitability in Fig. [2je), and signifies in- 
terference of the two decoupled eigenmodes on the detector. 
When driving the heptamer by a plane wave [Fig. |2jf), or- 
ange dashed line], we find that the same two modes cause the 
Fano-dip reported in literature 1 11 ]. The dark Fano-feature in 
extinction coincides spectrally with the asymmetric CL peaks, 
i.e. the asymmetric radiative LDOS enhancement, evident in 
local excitation. In contrast to local excitation, plane wave 
driving strongly drives the superradiant mode, causing the 
weakly excited subradiant eigenmode to appear as a narrow 
dip on a high background. We conclude that the measured 
redistribution of LDOS in CL data signifies the same interfer- 
ence mechanism as the Fano dip in extinction, though with 
very different superradiant and subradiant mode amplitudes. 

Far-field interference in eigen-basis. To gain further in- 
sight in the response of the individual modes, we focus on 
plane-wave driving and Fano-interference in the extinction 



cross-section (Text of the heptamer [black line Fig. [4ja)]. In 
Fig. |4|b,c) we plot the induced (complex) polarization p 
of the eigenmodes. Extinction (work done by the driving 
Im[Eeajt • p], normalized to incident intensity) can readily 
be split into contributions from different modes [Fig. |4ja)]. 
The total extinction cross-section is in excellent quantitative 
agreement with reported numerical results [11], underlining 
the suitability of a dipole model. The superradiant mode pro- 
vides a broad positive extinction [red line in Fig. |4ja)]. The 
Fano-dip in the sum is created by the subradiant mode [blue] . 
Its surprising negative contribution to dext indicates that the 
mode feeds energy back into the driving field, equivalent to 
destructive far-field interference. This energy cannot result 
from direct driving of the superradiant mode and subsequent 
amplitude transfer [27], since the eigenmodes are strictly de- 
coupled. Since the projection of plane wave driving on the 
eigenmodes has no frequency dependence, the excitation am- 
plitudes in Fig.|4jb) simply follow the eigenpolarizabilities in 
Fig.[3];c,d). On the Fano-dip, the subradiant mode's net dipole 
moment compensates that of the superradiant mode such that 
their superposition barely radiates. Importantly, the subradi- 
ant eigenmode must possess nonzero dipole moment (and ex- 
cellent radiation pattern overlap with the superradiant mode) 
to lead to a PIT dip in extinction. It is this nonzero dipole mo- 
ment together with a very strong, narrow eigen-polarizability 
[blue lines in Fig. [3];c,d)] that also ensures that the subradiant 
mode can give rise to an enhanced radiative LDOS at the cen- 
tral particle, as seen from Fig. |2jf). This eigenpolarizability 
exceeds the one of a single isolated particle and is narrower 
in width. We therefore draw the counterintuitive conclusion 
that introducing sub-radiant modes that are usually associated 
with dark PIT phenomena can actually enhance the capabil- 
ity of optical antennas to create bright and efficient emitters 
with a large radiative LDOS. The theoretical analysis tool that 
we have applied to the specific case of a heptamer will be of 
great value to explore the application of PIT to spontaneous 
emission enhancement. Our approach goes beyond both brute 
force numerical techniques and electrostatic hybridization and 
provides a generic, analytical framework that for the first time 
reveals the true eigenmodes of a structure exhibiting PIT, their 
polarizability, brightness, and response to any near- or far- 
field driving. 

Identification of bare and dressed states. Fano interfer- 
ence is inextricably linked to coherent transfer of amplitude to 
states that are not directly driven |27|. Indeed, this is the pic- 
ture in which the discussion of PIT has been led thus far 0|- 
[13 il . Necessarily, the eigenmode basis of any linear system 
never involves amplitude transfer, as eigenmodes are decou- 
pled. Figure |4ja-c) confirms this viewpoint, which is analo- 
gous to analyzing atomic FIT in terms of dressed states [IJ 
that decay to the same continuum with opposite phase. The 
alternative view on FIT, in which amplitude transfer does oc- 
cur, is as an interference of different pathways of coherently 
coupled bare states |1|. The basis of SALCs [Fig. [3^] pro- 
vides the analog of such bare states in FIT for the PIT sys- 
tem. The contributions of the three SALCs to (Text plotted in 
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Figure 4. (a) Calculated contribution of eigenmodes [Fig. [SJd)] to 
extinction cross section of plasmonic heptamers. (b) Eigenmode am- 
plitude and (c) phase relative to driving field. (d,e,f) Corresponding 
graphs in the non-diagonal bare state basis [Fig. [3^)]. 

Fig. [4jd) reveal a dip in the broad band of the super-radiant 
SALC [red line]. However, also the sub-radiant SALC [blue] 
contributes significantly. Fven though the contribution of the 
dark SALC [green] to (Text is strictly zero it is crucial for the 
Fano-dip. Figure |4je) shows that the dark SALC acquires an 
amplitude at the Fano-dip as large as that of the subradiant 
SALC. However, around the Fano dip the polarization p of 
the dark SALC is locked in amplitude and phase to the subra- 
diant SALC [Fig.|4];e,f)]. One can hence view the dark SALC 
as a screening effect that allows to describe the three- state in- 
teraction as the interaction between just two oscillators I 27 L 
with renormalized resonance frequencies and strengths. In 
conclusion, when taking the SALCs as a basis of bare states 
in which IM is not diagonal, amplitude transfer occurs to a lin- 
ear combination of the sub-radiant and dark SALC [Fig.|4je)]. 
The link between the complementary interpretations of FIT in 
scattering systems via either amplitude transfer or far-field in- 
terference is ultimately provided by the optical theorem which 
constrains IVE and thereby inextricably links radiated power 
and induced complex dipole moments to satisfy energy con- 
servation. 

Conclusions. In conclusion, we measured a spatial redis- 
tribution of the LDOS of plasmonic heptamers using CL at 
the spectral position that coincides with the reported Fano ex- 
tinction dip. Both LDOS redistribution and Fano dip result 
from interference in the far-field of the same two eigenmodes, 
excited in different coherent superpositions. Our findings 
bear a plethora of exciting prospects to harness the near field 
of plasmonic molecules, especially in the context of sponta- 
neous emission, or control of any process in sensing, spec- 
troscopy and non-linear optics that benefits from enhanced 
LDOS. Strikingly, optimizing a Fano dip in extinction re- 
quires engineering of radiation pattern overlap between two 
involved modes that ensures an enhanced LDOS by construc- 
tive interference. Also, completely dark SALCs may be uti- 
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lized to optimize resonances and interaction strengths, sim- 
ilar to the screening dark SALC in the heptamer. For flu- 
orescence appHcations, one typicaUy requires simultaneous 
optimization of pump field generated from a far field beam, 
and optimization of the LDOS at the redshifted emission fre- 
quency. Our generic quantitative solution approach may yield 
universal bounds on how to optimize the solution to such a 
problem using optical antennas with a Fano resonance. 
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